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Abstract 

In the Minimal Supersymmetric Standard Model, the masses of the charged Higgs boson (H^) 



■ and the CP-odd scalar {A) are related by M'^+ = M\+ml^ at the Born level. Because the coupling 

o ■ _ 

of W^-A-H^ is fixed by gauge interaction, the Born level production rate of qq' — > W^* — > AH^ 



depends only on one supersymmetry parameter - the mass (M^) of A. We examine the sensitivity 



^ ' of the LHC to this signal event in the j4(— > hb)H^{-^ t^v^) and A{-^ bh)H^{-^ th) decay channels. 

We illustrate how to test the mass relation between A and in case that the signal is found. If 
^ ' the signal is not found, the product of the decay branching ratios of A and predicted by the 



MSSM is bounded from above as a function of Ma- 
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I. INTRODUCTION 



One of the top priorities of current and future high-energy colhders, such as the Fermilab 
Tevatron and CERN Large Hadron Colhder (LHC), is to probe the mechanism of the elec- 
troweak symmetry breaking. In the Standard Model (SM) of particle physics, this amounts 
to searching for the yet-to-be-found Higgs boson. It is also possible that the mechanism of 
electroweak symmetry breaking originates from new physics beyond the SM. Very often, the 
low energy effective theory of new physics model predicts the existence of extended Higgs 
sector that contains more than one Higgs boson. Hence, in order to discriminate each new 
physics model from its alternatives, it is important to detect these additional Higgs bosons 
and to measure their properties. 

Supersymmetry (SUSY) is one of the most commonly studied new physics models. The 
Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) is known as a special 
case of the Two Higgs Doublet Model (THDM) with the type-II Yukawa interaction The 
THDM yields five physical scalar states, i.e., two CP-even [h nd H), a CP-odd {A) and a 
pair of charged Higgs bosons {H"^). 

In contrast to the free parameters in the Higgs sector of the THDM, the masses and 
mixing angles among the Higgs bosons in the MSSM are determined by the gauge couplings 
due to the requirement of supersymmetry. One of the striking features of the MSSM Higgs 
sector is that the mass {M^) of the lightest Higgs boson h is predicted to be smaller than 
the mass of the Z boson at the tree level and less than about 130 GeV after including 
contribution from radiative corrections Q]. The discovery of such a light Higgs boson could 
be a strong hint of the MSSM, as the mass of the SM Higgs boson has to be approximately 
between 130 GeV and 180 GeV for the SM to be a valid theory all the way up to the Planck 
scale (10^^ GeV) However, the discovery of a light Higgs boson does not exclude new 

physics models in which the light Higgs boson naturally has a mass less than 130 GeV, even 
when the Planck scale is taken as the cut-off scale of the models. Two examples are the 

n 

non-SUSY THDM and Zee model p. Therefore, a detailed study of the Higgs sector is 
necessary to discriminate the MSSM from other new physics models. Needless to say that a 
general test of the MSSM should also include the experimental identification of super-partner 
particles, such as sfermions, charginos and neutralinos. 
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Many studies have been done in the hterature to show how to detect heavier MSSM Higgs 
bosons produced at the LHC 0, Q]- A hght charged Higgs boson with mass Mj^± Knit—rrib 
may be produced via the top quark decay ^ or via the pair production qq H^H~ j^, 
gg H^H~ jl^ and qq H^H^qq 11 1. The main production channels for the heavier 



may be those associated with heavy quarks such as gh H^t 12 1 and qb q'hH^ llSll . 
The associated production of the W boson and charged Higgs boson 66, gg H^W^ |l4| 
may also be substantial, especially for small and large tan /3 values. The possible disco very 
channels from the decay of H are 

Furthermore, the CP-odd Higgs boson may be produced mainly in association with bottom 
quarks (qq.qq Abb) 1^. The possible discovery channels from the decay of A are A — >■ 

n n n n 

|2| , v4 — >• 77 and A — »• Zh [2i| for low tan (3 values, and A [2l| or rr [22| for 

large tan P values. The production and decay channels discussed above can be used to probe 
a large part of MSSM parameter space in the versus tan (3 plane, except for intermediate 
values of tan f3 and large values of M4 where only the lightest Higgs boson h is likely to be 
found. 

In a recent paper j^], two of us proposed to study the associated production of CP-odd 
scalar (A) and charged Higgs boson (H^) via qq' W^* — AH^ at hadron colliders to test 
the MSSM. It was shown that at the Born level the production of this process depends only 
on one SUSY parameter - the mass (M4) of A. This is because in the MSSM the masses of 
H"^ and A are related by the mass of the boson {mw) as^ 

M|± = Ml + ml, , (1) 

and the coupling of W~-A-H~^ is fixed by gauge interaction. In case that can be 
measured from data by reconstructing its decay particles (such as bb pair), the Born level 
production rate of this process is uniquely determined and does not depend on any other 
SUSY parameters. Therefore, studying this process can test the product of the Higgs boson 
decay branching ratios (in contrast to the product of decay branching ratios and production 
rate) predicted by the MSSM. The above conclusion holds for any SM decay model of A 



^ We note that in general, a CP violating phase can enter the Higgs sector of the MSSM, so that the 
CP-even Higgs bosons can mix with the CP-odd Higgs scalar and the mass relation Q does not hold any 
more. In this paper, we shall focus our study on the MSSM with a CP invariant Higgs sector. 
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and H"^. For example, if the decay mode of A — > 66 and r^Vr in the AH^ event 

is studied and no excess is found in experimental data for a given mass bin of (hence, 
Mfj±), one can constrain the MSSM by demanding the product of the decay branching ratios, 
Br(y4 —>■ bb) x Br(if+ r^Ur), to be bounded from above as a function of Ma- Similarly, 
applying the same strategy, one can constrain the product Br(y4 —>■ X) x Bt{H^ —>■ Y) 
for any decay modes X and Y predicted by the MSSM, as a function of only one SUSY 
parameter - Ma- 

Another interesting feature of this process is that the kinematic acceptance (therefore, 
the detection efficiency) of the signal events does not depend on the choice of the other 
SUSY parameters because both A and are spin-0 (pseudo-) scalar particles so that the 
kinematic distributions of their decay particles can be accurately modelled. (It is an isotropic 
distribution in the rest frame of A or H^-) Moreover, the one-loop electroweak corrections 
to the mass relation and the W~AH~^ coupling are generally smaller than the other 
theoretical errors (such as the parton distribution function uncertainties) and the expected 
experimental errors (such as the mass resolution of Higgs boson decaying into jets). Despite 
that for some extreme choice of the MSSM parameters, the mass relation ffl can receive a 

, its effect to the AH'^ production rate is found to be small. 
Therefore, in this work we shall only perform a Born level Monte Carlo analysis to study 
the sensitivity of the LHC to this production channel. 

To detect the signal event at hadron colliders, it is necessary to suppress the potentially 
large background rates predicted by the SM. In this paper, we show how to detect such a 
signal at the LHC for the specific decay mode of AH"^ bbr^Ur with r+ tt+z/^. While it 
is possible to include other decay channels of AH"^, we shall restrict ourselves to this simple 
decay mode and use it as an example to illustrate the idea of Ref. I23I. Namely, studying 
the associated production of A and at hadron colliders can discriminate the MSSM from 
its alternatives, e.g., a general two-Higgs-doublet model (THDM), and test the prediction 
of the MSSM on the product of the decay branching ratios of A and H^- 

The rest of the paper is organized as follows. In Sec. II, we briefiy review the result of 
Ref. for the production rate of the signal event at the Fermilab Tevatron and CERN 
LHC up to the one-loop order QCD and electroweak corrections. In Sec. HI, we present a 
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FIG. 1: The LO (dotted lines) and NLO QCD (solid lines) cross sections of the AH+ and AH' 
pairs as a function of at the Tevatron (a 1.96 TeV pp collider), and the LHC (a 14TeV pp 
collider). The cross sections for AH^ and AH^ pair productions coincide at the Tevatron for 
being a pp collider. 

Monte Carlo analysis for studying the sensitivity of the LHC to the detection of the signal 
event. We show that all the dominant SM backgrounds can be effectively suppressed by 
making use of the polarization of r^. Sec. IV contains our conclusions. 

II. PRODUCTION CROSS SECTION OF SIGNAL 

At hadron colliders, the Born level process for the production of a AH^ pair is qq' — 
W^* AH^. The coupling of W^-A-H^ is induced from the gauge invariant kinetic term 
of the Higgs sector [l| : 



and its strength is determined by the weak gauge coupling g. Hence, the Born level cross 
section a{pp, pp AH"^) depends only on the masses of A and if^. Its squared amphtude, 
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(2) 



■int 
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after averaging over the spins and colors is 

3 (s2 -m^)2 + m^r^ 

where P = - M\ with = (s + M| - Mjj+)/{2^) and ^ is the polar angle of 

A in the center-of-mass (cm. ) frame of A and H^. In the MSSM, the masses of the 
charged Higgs boson H"^ and the CP-odd scalar A are strongly correlated. At the Born 



level, Mh± = a/ M\ + m^. In Fig. ^ we show the inclusive production rate of AH^ as a 
function of for the Tevatron (a 1.96 TeV proton-antiproton collider) and the LHC (a 
.4TeV proton-proton collider). Here, the CTEQ6M parton distribution functions (PDF) 
2^ are used and both the renormalization and the factorization scales are chosen to be 
the invariant mass (v^) of the AH"^ pair. We note that for the ch — > AH^ subprocess, 
in addition to the Cabbibo-Kobayashi-Maskawa suppressed s-channel VT-boson diagram, 
there is a t-channel diagram whose contribution becomes large for a large value of tan j3. 
Nevertheless, the cb AH^ contribution to the inclusive AH^ rate is small. For example, 
its contribution to the total rate is less than 0.01% and 0.1% at the Tevatron and the LHC, 
respectively, for tan /? = 40 and = 90 GeV. For a smaller value of tan /3, its contribution 
becomes negligible. Hence, we shall ignore its contribution in the following discussion. 

It is easy to include the next-to-leading order (NLO) QCD correction to calculate the 
production rate of the signal event because at the one-loop order only the initial state 
receives radiative contribution which is similar to the NLO QCD correction to the IV-boson 
production, but with a different invariant mass, at hadron colliders. The result is also shown 
in Fig. H We find that typically the NLO QCD production rate is about 20% higher than 
the leading-order (LO) rate when the same PDF is used. The higher order (af or above) 
QCD correction is estimated to be about 10% at the Tevatron and less than a percent at the 
LHC, for = 120 GeV, when the factorization scale is varied around the AH^ invariant 
mass by a factor of 2. 

The dominant one-loop electroweak corrections to the qq' —>■ AH^ process come from the 
loops of top (t) and bottom {b) quarks as well as their supersymmetric partners, i.e. stops 
(^1^2) and sbottoms (61,2), in the MSSM, due to their potentially large couplings to Higgs 
bosons. 

In the following, we discuss the quark- and squark-loop radiative corrections to the ef- 
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fective coupling of W^AH^ and to the mass relation (^J. In our calculation, we adopt the 
on-shell renormalization scheme developed by Dabelstein in Ref. 3|. (For completeness, 
we outline the renormalization scheme in Appendix A). 

The part of the one-loop effective coupling of W^AH^, that is relevant to the production 
process qq' —>■ H^A, can be written as^ 

Mi^HAiQ') = -f(P^-P^r [l + F^'Hq')] , (4) 

where g^, and p'^ are the incoming momenta of W~^, A and H~ , respectively, and g is 
the effective weak gauge coupling evaluated at g^. Hence, the radiative correction to the 
cross section of the sub-process qq' —>■ AH^ at the one-loop order is 

= 2ReF«(g2), (5) 

where the detailed calculation for F^^^q"^) is summarized in Appendix B. As shown in 
Eqs. (I2HI), and (jKTj) . the quark-loop contribution is proportional to the squared Yukawa 
coupling constants yti= 2mf cot"^ f^/v"^) and y^(= 2m^ tan^ ^). In the large rrit or large 
m;,tan/3 limit, it can be written as 



quark IQjf'^ 



(6) 



where Nc{= 3) is the number of colors. Since and y'^ are at most 0{1) for tan/3 ^ 1 and 
mt/mb, respectively, -Fq^ark at most a few percent for 1 ~ tan/5 ~ nit/mb. 

We also calculate the squark-loop contribution. As compared to the quark-loop effects, 
the squark-loop effects are rather complex because of the additional SUSY parameter de- 
pendence. The mass eigenstates /i^2 {f = t or b) of the squarks are obtained from the weak 
eigenstates /l^r by diagonalizing the mass matrices defined through 
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where, M| = M?+m^ + (m| cos2/5)(T/^ -Q/s^) and M| = M? ^+m^ + (m| cos2/3)Q/s^. 
In this expression, M? M?. (for f = t) and M"-^, (for f = b) are the soft-breaking masses for 



The other form factor, (p^ does not contribute to this process for massless quarks. 
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fi, iji and respectively; sw = sinOw with 9\y being the weak mixing angle; Tf^ is the 
isospin of the left-handed quark fi] and Q/ is the electric charge of the quark /. Moreover, 
Xf = At — yUcot/9 and Xf, = Af, — fitanf], where At (Ah) is the trilinear A-term for t (b), 
and fi is the SUSY invariant higgsino mass parameter ^27]. For completeness, we have listed 
all the relevant squark and Higgs boson couplings in Appendix C, so that the squark-loop 
contributions to F^^\s), cf. Eq. can be directly read out from the results in Eqs. (j3U|) . 
dSHl) and (f^. 

To examine the effect of one-loop electroweak corrections, we shall discuss two limiting 
cases below. Firstly, we consider the cases with fi = At = Af, = 0, i.e., the cases without 
stop mixing {\Xt\ = 0) and sbottom mixing {\Xh\ = 0). Under this scenario, the masses 
of squarks are proportional to the typical scale of the soft-breaking mass, denoted as M. 
Since all the relevant couplings between squarks and Higgs bosons are independent of the 
soft-breaking masses Mq, and (see Appendix C), the squark-loop effect is decoupled 
and its contribution is very small for a large value of M, where M = Mq ~ Mq- ~ M^. 
For a smaller M, -^squark becomes larger. However, M cannot be too small because a small 
M implies light squarks whose masses are already bounded from below by the direct search 
results 13]. Furthermore, as to be shown later, the case with a small M is also strongly 
constrained by the p parameter measurement. Secondly, we examine the case with a large 
stop mixing, assuming mt|Xj| ~ ^ m|. Such a large stop mixing leads to a large mass 
splitting between ii and ^2 so that m^^ ~ 0{mz) and m^^ — V^M, while m^^ ^ ~ M. The 
leading squark contribution to F^^\q'^) can be expressed as 



F 



(1) _ , -N, 



squark 



167r2 



4 \M \ 6 \M 



(8) 



with Yi = ^ (At cot (3 + /i) and = ^{Abtan(3 + /i). Since in this case |A| ^ \M'^/mt ± 
cot/?|, we have \Yi\ ~ 0{M^/v) for | /i| ~ M and tan/5 ~ 1. When \Ah\ ~ \At\ > \p\ and 
tan/5 ~ rrit/mj,, we find \Yi\ ~ 0{M^/v). Thus, with a large stop mixing (m(|Xt| ~ M^), 
-^squark proportioual to the soft-breaking mass scale M, and does not decouple in the large 
M limit. However, such kind of model is strongly constrained by the p-parameter (or the 
T-parameter) data. With a large stop mixing (M^ ~ mt\Xt\), the squark contribution to 
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the p-parameter (cf. Appendix D) is 



Apsquark^ (2.2 X 10-3) — . (9) 



Since any new physics contribution to the p-parameter has to be bounded by data as 

- 1.7 < Apnow X 10^ < 2.7, at 2a level, (10) 

the scale M cannot be too large in this case. Consequently, the above -^squark constrained 
to be smaller than a few percent as long as p^ is not much larger than M^. 

To examine the effect from the stop and sbottom loops to the production rate of AH^, 
we consider 4 sets of SUSY parameters, as listed in Table HI which give the largest allowed 
deviation in the p-parameter. Set 1 and Set 2 represent the cases without either a stop mixing 
{Xt = 0) or a sbottom mixing {Xt = 0), and Set 3 and Set 4 are the cases with a large 
stop mixing (mi|Xf| ~ M^) and m^^ ^ 100 GeV. The K^^\s) factor, as defined in Eq. 
is shown in Fig. |21 as a function of the invariant mass (^/s) of the constituent process for 
= 90 GeV. It is clear that the quark-loop contribution to K^^\s) dominates the squark- 
loop contribution. For these four sets of SUSY parameters, the squark-loop contribution is 
smaller than the quark-loop contribution by about a factor of 100. Generally, the squark 
contributions are at most a few percent, unless |p| is taken to be very large as compared to 
the scale M. We have checked that this conclusion does not change when the assumption 
Mq ~ M^ ~ M|, is relaxed to some extent. Including both the quark- and squark-loop 
contributions to K^^\s), we found that the correction to the hadronic cross section of H^A 
production in the invariant mass region just above the H~^A threshold, where the constituent 
cross section is the largest, is at a percent level. In summary, to be consistent with the low- 
energy data and the direct search results for stops and sbottoms, the one-loop electroweak 
correction to the production rate of pp, pp AH^ is small (at most a few percent). 

Next, we discuss the one-loop corrections to the mass relation (^. Let us parameterize 
the deviation from the tree-level relation by 6, so that at the one-loop order 



Mh± = sjMl + ml, {1 + 5). (11) 

We note that in our renormalization scheme (see Appendix A), and my/ are the input 
parameters, but Mh+ is not. The one- loop corrected mass of the charged Higgs boson M^± 



TABLE I: The SUSY (input and output) parameters used in Fig. |21 





Setl Set2 


Set 3 


Setl 


Mq = M^ = (GeV) 


106 


84 


408 


409 


tan P 


2 


40 


2 


40 


At = Ab (GeV) 








+1261 +1119 










+300 


+300 


(GeV) 


197 


184 


100 


100 


^i, (GeV) 


199 


188 


612 


611 


m-,^ (GeV) 


108 


88 


407 


373 


m-,^ (GeV) 


116 


103 


412 


447 


Apsquarks X 10^ 


2.72 


2.70 


2.71 


2.70 



can be obtained by solving 



= Det 



(2) 

G+G~ 
(2) 

H+G- 



(p^) r 



(2) 

G+H^ 
(2) 



(12) 



where T^^^ (p^) represent the renormahzed two-point functions in the basis of the renormahzed 
Goldstone boson (G^) and charged Higgs boson {H"^) fields. Here, the notation "Det" 
denotes taking the determinant of the 2x2 matrix. One of the solutions of the above equation 
is p"^ = 0, which corresponds to the charged Nambu-Goldstone mode, and another is M^±. 
At the one-loop level, the pole mass of the charged-Higgs boson, in our renormalization 
scheme, can be calculated from 



M^± = Ml + 

+UAAiMl) - Uh+hAMa + mli) + nw.iy(m^), (13) 

where n<^<^(g2) (0 = A,i7^,and W) are the self-energies. For completeness, we summarize 
the quark and squark contributions to the self-energies of A and in Appendix E. 

When y4(.ft and yU are zero (i.e., no- mixing case), the leading contribution to 5 is found to 
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FIG. 2: The K-factoT, K^'^\s), of the constituent process qq' H+A for = 90 GeV, as a 
function of the invariant mass ^/s of qq'. The soHd hues correspond to the top and bottom quark 
contribution. The cases where the squark-foop contribution is included are described by the dotted 
hues for those without stop mixing (Set 1 and Set 2) and by dashed hnes for those with maximal 
stop mixing (Set 3 and Set 4), respectively. 



be 



6 



2 2 



In 



M2 



(14) 



87r2t;2 \Ml + ml,) sin^ /5 cos2 /5 

This correction is substantial for tan/? ~ rrit/mh and ^ mf. Applying Eq. p3|) with 
the complete expression of U^^^q"^), we found 6 to be less than 4.9% for 2 < tan/3 < 40, 
M < 2 Te V and Ma > 90 GeV. Our result agrees well with Ref . [s^ , in which an approximate 
formula was presented for ^ mf. 

For the cases with nonzero Ah t and fi parameters, 6 receives extra contributions, which 
are proportional to A^f^/M"^, A^^ii^ /M'^ and /i^/M^ 30, |^| originated from the squark 
couplings [cf. Eqs. (jHUj) . (jH7|) . (jlOj) and PTj) ] and squark masses. For the Set 3 and Set 4 
parameters listed in Table HI 5 is found to be less than 5.3% and 3.6% for Ma > 90 GeV, 
respectively. In summary, as long as |A,f)| and |yu| are not too large as compared to M, 
in a wide range of the parameter space that is allowed by the available experimental and 
theoretical constraints, 6 does not exceed 10%. (In some corner of SUSY parameter space, S 
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can be larger than 10% [2^. Nevertheless, the AH~^ production rate at the LHC will not be 
strongly modified as long as is not too large so that the mass threshold effect becomes 
important.) 

Supported by the finding that the one-loop electroweak corrections to the W^AH^ 
coupling and to the mass relation = M\ + are generally smaller than the 

other theoretical errors (such as the parton distribution function uncertainties) and the 
expected experimental errors (such as the mass resolution of Higgs boson decaying into jets, 
which is about 10 GeV for a 100 GeV Higgs boson), we conclude that studying the process 
pp ,pp ^ W^* AH^ allows us to distinguish the MSSM from its alternatives by verifying 
the mass relation (^, or Eq. (fTT|) . when the signal is found. If the signal is not found, 
studying this process can provide an upper bound on the product of the decay branching 
ratios of A and as a function of only one SUSY parameter - M^. 



III. MONTE CARLO ANALYSIS 



In this section we discuss how to detect such a signal at hadron colliders. As a concrete 
example, we examine the specific decay channel of AH^ hhr^Vr with r+ tt+P^. As 
shown in Fig. ^ the signal cross section at the Tevatron is much smaller than that at the 
LHC. Hence, we shall first discuss its phenomenology at the LHC, and then comment on the 
result at the Tevatron. Furthermore, since the NLO QCD and electroweak corrections to 
the signal cross section are small, for simplicity, we shall perform the Monte Carlo analysis 
at the Born level. 

In the MSSM, the mass of the heavier CP-even Higgs boson [H) is about the same 
(approximately within 10 GeV) as Ma for Ma ~ 120 GeV and tan/3 ~ 10. In this case, 
gg' HH^ can produce similar final states as qq' AH^. In contrast to the coupling 
of W^AH"^ whose strength is ^, the tree level coupling of W^HH"^ has the strength of 
— sin(/3 — a), though with the same momentum dependence. Since sin^(/5 — a) ~ 1 for 
Ma ~ 190 GeV and tan/5 ~ 10, the production rate of HH^ is almost the same as AH^ in 
the decoupling limit. Thus, for a large Ma, when both of them decay into the same decay 
channels, it will be very difficult to separate the production of AH"^ from HH^ unless a 
ffne mass resolution can be achieved experimentally. Nevertheless, studying different decay 
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channels can help to separate these two production modes. For example, a heavy H can 
decay into a ZZ pair at the Born level, but A cannot. In the following analysis, we shall 
consider both the AH^ and HH^ production channels that contribute to the same final 
state. 

Either in the MSSM or the Type-11 THDM, for a large tan (3 value, the dominant decay 
mode is h,H,A—* bb for the neutral Higgs bosons {h, H, A), and r+z/ for a charged 

Higgs boson whose mass is smaller than the top quark mass. Due to the missing energy 
carried away by the final state neutrino, it is not possible to directly reconstruct the mass 
of if"*" in the decay mode t^v. But, the transverse mass of the charged Hi ggs boson 

can be reconstructed from the r jet momentum and missing transverse energy [l^. For 
mH± ^ 200 GeV, the dominant decay mode is tb, in which Mh+ can be reconstructed 

after properly choosing the longitudinal momentum of the neutrino (with a two-fold solution) 
from t decay. To cover both decay modes in our study, we consider the following three 
benchmark cases with tan/? = 40: (i) Ma = 101 GeV (and Mfj+ < rrit + mh), and tu 
being the dominant decay mode; {ii) Ma = 166 GeV (and Mh+ ~ + m^), and if + tu 
being the dominant decay mode; [iii) Ma = 250 GeV and tb being the dominant 

decay mode. 



A. Ma = 101 GeV 

Taking Ma = 101 GeV, tan/3 = 40, and all the other SUSY parameters (soft-breaking 
masses and /i-parameter) to be 500 GeV, one can calculate the masses, widths and decay 
branching ratios of the Higgs bosons using the numerical program HDECAY 0- We 
find that Mh = 97 GeV, Mh = 113 GeV and Mh± = 126 GeV, and the relevant decay 
branching ratios are B{A -> bb) = 0.91, B{H bb) = 0.90, B{H+ r+z/) = 0.98, and 
B{t^ TT^u) = 0.11. The relevant total decay widths are = 3.7 GeV, Th = 0.38 GeV 
and Th± = 0.43 GeV. For clarity, we summarize the above results in Table IHI In spite of 
the small Mh, this set of SUSY parameters is compatible with CERN LEP data in searching 
for neutral Higgs bosons of the MSSM js^, and yields sin(/? — a) = 0.33. 



13 



TABLE II: The mass, width, production rate and decay branching ratio of Higgs bosons used in 
the Monte Carlo analysis. Here, tan /3 = 40, and all the other SUSY parameters (soft-breaking 
masses and //-parameter) are taken to be 500 GeV. 



Sets 


A 


B 


C 


niAjVA (GeV) 


101 / 3.7 


166 / 5.6 


250 / 7.9 


mhl^h (GcV) 


97 / 3.3 


112 / 0.04 


112 / 0.01 


mnlVn (GeV) 


113 / 0.38 


163 / 5.5 


248 / 7.8 


m,^+/r^+ (GeV) 


126 / 0.43 


182 / 0.68 


261 / 4.2 


a^AH^) fb 


164 


36.5 


5.4 


a{HH+) fb 


15.0 


37.2 


5.4 


B{A-^ bb) 


0.91 


0.90 


0.89 


B{H ^ bb) 


0.90 


0.90 


0.89 


B{H+ T+v) 


0.98 


0.90 


0.21 


B{t+ TT+U) 


0.11 


0.11 


0.11 


B{H+ tb) 


0.00 


0.09 


0.79 


B{t^ be+v) 


0.11 


0.11 


0.11 


sin(/3 — a) 


0.33 


0.997 


0.999 



1. Signal and background Processes 

In this case, the signal event is produced via 

g + g' ^ A{^hh)H+{-^vT+{-^T:^V)), 

q + q' ^ H{-^ bb)H+{^ z/r+(^ 7r+z>)) . (15) 

The experimental signature of the signal events is the production of ir'^bb associated with 
a large missing transverse momentum carried away by neutrino and anti-neutrino. In this 
study we shall assume that both b and b jets in the signal event can be tagged with a total 
efficiency of 50%, so that the four dominant SM background processes are: 

Wbb : qq bbW^, W+ t+v, r+ ^ vr+i/ , 
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tb : qq' ^W* ^tb,t^ W+ t+v, t+ ^ tt+P , 

Wg: qg^ q'tb, t W+b, W+ r+u, t+ ^ tt+P , 
ti: qq, gg^tt,t^W+b,W+ ^r+u,T+ ^7r+iy,t^bW- . (16) 

The first two processes are the intrinsic background processes generated by the SM, and 
yield the same experimental signature as the signal process at the parton level. In addition 
to the Wbb and tb processes, the other two important SM background processes are Wg and 
tt. When the forward jet {q') in the Wg fusion process is not detected (either faUing into a 
large rapidity region or carrying a too small transverse momentum to be detected), it will 
mimic the experimental signature of the signal event. Similarly, when the decay products of 
W~ in t — > bW~ escape detection, the tt event yields the same event signature as the signal. 
To study how often this situation happens, we separately consider the leptonic and hadronic 
decay modes of W~ at the parton level. It is straightforward to study the decay modes 
W~ — > £~Pe with £~ — e' and For W~ — > t~Vt, we consider the decay of t~ into the 
£~P^Pt-, 7r~P^ or p~Vj- mode. We find that the dominant contribution after imposing the veto 
cuts (to veto events with extra charged lepton with large transverse momentum or hadronic 
activities in the central rapidity region) comes from W~ i~Pi and W~ p^r^(-^ f^z/^P^), 
with i~ = e~ and Since in the end of analysis, the ti rate is smaller than the Wbb rate 
by more than one order of magnitude, we shall defer the detailed discussion on this part of 
the analysis to Sec. III.C, where Case C is considered. As to be shown later, for Case C, the 
dominant background comes from the tibb events which mimic the signal event signature 
when the decay products of the extra W~, from i decay, escape detection. 



2. Analysis 

a. Basic cuts: To compare the relevant background event rates to the signal event 
rate, we shall assume the integrated luminosity of the LHC to be 100 fb~^, and require the 
transverse momentum (p^) and the rapidity (77^) of b, b and tt"*" to satisfy the following basic 
cuts: 

> 15 GeV, \r]'^\ < 3.5, AR > 0.4 . (17) 
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Here, we have assumed a perfect detector that can precisely measure the four-momenta of 
the final state partons and require the separation in AR (= a/ (A0)^ + (Ar^)^) between any 
two observable final state partons (not including neutrinos) to be larger than 0.4, where Ac/) 
and Arj are the separation in azimuthal angle and rapidity, respectively. (We shall comment 
on the effect due to the finite resolution of the detector in the later part of this section.) 
Furthermore, in the Wg event, the additional q' jet (preferably in the forward direction) is 
required to escape detection, i.e., either its transverse momentum is less than lOGeV or its 
rapidity (in magnitude) is larger than 3.5. Similarly, the charged lepton i~ (= e~,fi~,r~) 
from the decay of i in the ti event is required to be undetected, i.e., < 10 GeV or 
|?7^ I > 3.0. For clarity, we shall only include the positively charged state (i.e., vr"'') in the 
following discussion. The inclusion of the opposite charged state (i.e., vr") and other decay 
modes of will be discussed in Sec. Ill D. The numbers of the signal and background 
events after imposing the above cuts are summarized in the second column of Table IIIIl in 
which the 6-tagging efficiency (50%, for tagging both b and b jets) is included. The last 
three rows in Table lllll show the ratio of signal (S) to background (B) event number, the 
statistical significance of the signal, and the statistical uncertainty in the measured signal 
event rate. 

At this stage of the analysis, the background rate is an order of magnitude larger than the 
signal rate. Moreover, the dominant background comes from the Wbb process, followed by 
the s-channel single-top {tb) process, the t-channel single-top {Wg) process and the tt pair 
process. To detect the signal event, further kinematic cuts are needed. In order to study 
the efficient cuts that can largely suppress the background rates while keeping most of the 
signal rates, we examine the distributions of the missing transverse energy (-^t), transverse 
momentum of b jet (p^), transverse momentum of vr"*" (pj-), and the invariant mass of bb 
pair {fTij^). Their distributions are shown in Fig. El (The vertical axis shows the number of 
events per bin for all the figures presented in this section.) 

h. Missing As shown in Fig. |21 the typical size of the missing transverse energy 

in the signal event is larger than that in the Wbb, tb and Wg background events. This is 
because the main source of the missing energy in the signal and background events comes 
from the neutrino in the decay of and W^ , respectively, and the mass of considered 
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TABLE III: Numbers of AH'^ signal and background events for Case A, with Ma = 101 GeV in the 
bbir^^T channel, at the LHC with an integrated luminosity of 100 fb^^. The 6-tagging efficiency 
(50%, for tagging both b and b jets) is included, and the kinematic cuts listed in each column are 
applied sequentially. 





Basic Cuts 


> 50GeV 


> 40 GeV 


90 < A% < 110 GeV 


With smearing 


AH+ 


507 


391 


241 


216 


202 


HH+ 


48 


38 


24 





4 


Wbb 


11555 


3111 


864 


67 


62 


tb 


1228 


614 


163 


12 


11 


Wg 


567 


236 


68 


11 


9 


tt 


110 


80 


17 


2 


2 


Signal (5) 


507 


391 


241 


216 


202 


Bckg {B) 


13507 


4078 


1135 


92 


87 


S/B 


0.038 


0.095 


0.212 


2.35 


2.32 


s/Vb 


4.36 


6.12 


7.14 


22.5 


21.6 


VS + B/S 


0.23 


0.17 


0.15 


0.08 


0.08 



here is larger than the mass of the boson. The average value of fir in the ti event 
is larger than the other background processes, for the undetected charged lepton i~ from 
the decay of i also contributes to the missing energy of the event. The same argument also 
explains why the average missing energy in the Wg event (with an extra q' parton jet) is 
slightly larger than that in the Wbb and tb events. In Table IIIIl we show the number of 
signal and background events after demanding the missing energy of the event to be larger 
than 50 GeV. This cut increases the signal-to-background ratio by about a factor of 2.5 while 
keeping about 77% of the signal rate. The biggest reduction in the background rate comes 
from the Wbb event (reduced by a factor of 3.7). 

c. Transverse momentum of the ir-jet: Another big difference between the signal and 
the background event signature is in the distribution of the transverse momentum of vr"*" 
originated from the decay of r"*", cf. Fig. El It can be understood as follows. In the 
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TABLE IV: Numbers of HH^ signal and background events for Case A, with Mh = 113 GeV in 
the bh-K^I^T channel, at the LHC with an integrated luminosity of 100 fb^^. The 6-tagging efficiency 
(50%, for tagging both h and h jets) is included, and the kinematic cuts listed in each column are 
applied sequentially. 





Basic Cuts 


> 50GeV 


p'^ > 40 GeV 


105 < M^^ < 125 GeV 


With smearing 


HH+ 


48 


38 


24 


24 


22 


AH+ 


507 


391 


241 


26 


43 


Wbb 


11555 


3111 


864 


58 


54 


tb 


1228 


614 


163 


11 


11 


Wg 


567 


236 


68 


13 


12 


ti 


110 


80 


17 


2 


2 


Signal (5) 


48 


38 


24 


24 


22 


Bckg {B) 


13966 


4431 


1352 


110 


120 


S/B 


0.003 


0.008 


0.018 


0.22 


0.19 


s/Vb 


0.41 


0.57 


0.65 


2.26 


2.03 


VS + B/S 


2.47 


1.76 


1.55 


0.49 


0.53 



signal event, the coupling of H^r^Ur is almost left-handed due to the enhancement from 
large tan/3. Consequently, the r+ from decay is produced preferentially left-handedly 
polarized, cf. Fig. |3 On the contrary, in the background events, comes from the decay 
of the weak gauge boson via left-handed weak current. In the dominant background 
{Whh) event, the boson is predominantly left-handed polarized. Therefore, r+, as an 
anti-fermion, is right-handedly polarized when ignoring the small mass of r+. In the weak 
decay of r"*" — > tt^Vt-, the left-handed weak current forces Vr to be purely right-handed in the 
limit that the neutrino mass is ignored. Hence, to preserve angular momentum conservation, 
7r+ would prefer moving along the direction of a left-handed r"*", and against the direction 
of a right-handed r"*", cf. Fig. |3] Namely, the angular distribution of vr^ is (1 + cos6',r) in 
the rest frame of a left-handed r"*", and (1 — cos^^^) for a right-handed, where 6*,^ is the polar 
angle of vr"*" momentum with respect to the moving direction of r"*" in the rest frame of 
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FIG. 3: Distributions of p\<, and m^^ at the LHC for Case A, with Ma = 101 GeV and 
Mh = 113 GeV in the bbTr~^^T channel, after imposing the basic cuts specified in Eq. H17() . 

(or W^). Hence, the transverse momentum of vr"*" in H'^ z/^r+(^ tt+z/t-) is typically larger 
than that in — > UrT^i^ vr^z/^). To further suppress the background rate, we require 
> 40 GeV. This cut increases the signal-to-background ratio by a factor of 2.2 at the cost 
of the reduction in signal rate by 40%, while the background rate is reduced by a factor 
of 3.6. As clearly indicated in Table 11111 a similar reduction factor applies to all different 
background processes. 

d. Invariant mass of the bb-jet pair: At this stage of the analysis, the dominant back- 
ground rate still comes from the Wbb process whose rate is about 3.6 times the AH^ signal 
rate. Since in the Wbb process, the bb pair originates from a virtual gluon conversion, the 
invariant mass of bb pair {nif^i) is generally small. This is clearly illustrated in Fig. El A 
similar, though less dramatic, result after imposing the kinematic cuts > 50 GeV and 
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FIG. 4: r"*" is left-handedly polarized in t~^v, and right-handedly polarized in — > t^v. 

Moreover, 7r+ momentum depends on the polarization of t+. (The thin arrow represents the 
moving direction, and the bold arrow represents the spin direction of the particle.) Similar plots 
for are also shown in the right four diagrams. 

> 40 GeV, is also shown in Fig. for comparison. To further improve the signal-to- 
background ratio, we require the invariant mass of hh pair to be within ± a, where 
o" = 10 GeV is the expected experimental resolution for a 100 GeV Higgs boson decay- 
ing into a hb pair [s^ . As shown in Table IIIH the Whh background rate is reduced by 
a factor of 13, and the other background rates are also reduced by an order of mag- 
nitude. This yields S/B = 2.35 and S/\fB = 22.5, with a total of 216 signal events, 
q + q' ^ A{—* bb)H^{^ ut^{^ tt+z/)) , produced at the LHC with an integrated luminosity 
of 100fb~\ 

Up to now, we have focused our discussion on the AH~^ signal channel. A similar analysis 
can also be performed for the HH^ signal channel, whose result is shown in Table HVl where 
we have required the invariant mass of hh pair to be within Mh ± a with a = 10 GeV. As 
shown in Table UTl since the coupling of W^HH^ is suppressed by a factor of sin(/3 — a) as 
compared to that of W^AH^ and Mh > Ma, the production rate of HH^ is expected to 
be smaller than the AH^ rate by more than a factor of sin^(/5 — a). For M^ = 101 GeV and 
tan/3 = 40, this reduction rate is about 11. Consequently, this HH^ signal is difficult to 
detect, for the background rate is still larger by a factor of 5 even after requiring > 50 GeV 
and pl^ > 40 GeV, cf. Table HVl and Fig. El 
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FIG. 5: Distributions of m^^ for Case A, with Ma = 101 GeV and Mh = 113 GeV in the hhii+^T 
channel, after imposing the additional cuts: I^t > 50 GeV and > 40 GeV. 

We therefore conclude that for a small Ma, the mass of A can be determined from the 
Mf,5 distribution, as shown in Fig. El However, to measure the mass of H via this channel is 
challenging. 

e. Transverse mass of the charged Higgs Boson: As discussed above, the mass of A 
can be determined from the invariant mass distribution of the bb pair. In order to test the 
mass relation, Eq. (0), predicted by the MSSM, the mass of should also be measured. 
In the signal event, the decay product of the charged Higgs boson contains the charged 
pion TT"*" and the two (anti-) neutrinos (z/,- and Ur), which contribute to the missing transverse 
energy Hence, we can examine the transverse mass of defined as 

mr = '^2pJ^t(1 - cos0) , (18) 

where is the azimuthal angle between the 7r+-jet and the missing transverse energy It 
is expected that the shape of the ttit distribution should exhibit a Jacobian peak behavior 
when the transverse momentum of vr"*" is large, for the dominant missing energy comes from 
Ut- in u^r^i^ n^Ur)- This feature is clearly illustrated in Fig. IHl The signal event 

shows a Jacobian peak around Mh+ (126 GeV), while the dominant Wbb background event 
shows a Jacobian peak around mw (80.2 GeV) but with a long tail into the larger M^g region. 
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FIG. 6: Distributions of the transverse mass my, defined in Eq. H18() . for Case A, after imposing 
the additional cut: \mf,i — ""T-A/i^l < 10 GeV. 

The long tail comes from the Wbb events in which the transverse momentum of the 
boson is large. Since the other background event rates are much smaller than the AH~^ 
signal rate, we conclude that the charged Higgs boson mass can be accurately extracted 
from the transverse mass distribution. 

/. Detector Effects: In reality, the performance of the detector is not perfect. To study 
the effect due to the finite detection efficiency of the detector, we repeat the above Monte 
Carlo analysis after smearing all the final state parton momenta by a Gaussian distribution 
with 

AE _ 50% 

where E is the energy of the observed parton and the resolution of the energy measurement 
is assumed to be 50%\/E. The distributions of invariant mass m^j^ and the transverse mass 
rriT become slightly broader, as shown in Fig. [7| However, both the signal and background 
rates, as shown in the last column of Tables UTTl and HVl are almost the same as those obtained 
with a perfect detector. 

Finally, we note that in the above analysis we have limited ourselves to a very minimal 
set of kinematic cuts to enhance the signal-to-background ratio (by about a factor of 60) 
while keeping most of the signal events (about half of them). There are other distributions 
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FIG. 7: Distributions of invariant mass m^-^ (after imposing the basic cuts, and pl^ cuts) and 
transverse mass rriT (after imposing all the kinematic cuts) for Case A, with Mj^ = 101 GeV and 
Mh = 113 GeV, after smearing the observable parton momenta by 50% V^- 
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FIG. 8: Transverse momentum distribution of the bb pair for Case A, after imposing the additional 
cut: \m^i — rriyi/Hl < 10 GeV. 



that can further distinguish signal from background events. For example, as shown in Fig. 
the transverse momentum of the bb pair in the AH~^ signal events is typically larger than 
that in the background events, even after imposing the kinematic cuts > 50 GeV and 
> 40 GeV. This is because the scattering amplitude of the AH^ signal event is a p-wave 
amplitude, cf. Eq. (3), and Mh+ is larger than M^, so that A is boosted to produce a large 
transverse momentum which is typically about 100 GeV in this case. 
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B. Ma = 166 GeV 



We shall extend our analysis to the production of a charged Higgs boson with its mass 
only slightly more than the top quark mass (174 GeV) so that the dominant decay mode of 
remains to be the r+z/ channel. We take M4 = 166 GeV, tan (3 = 40, and all the other 
SUSY parameters (soft-breaking masses and /i parameter) to be 500 GeV. Consequently, 
the masses of the other Higgs bosons are calculated from HDECAY [32j , which gives M/j = 
112 GeV, Mh = 163 GeV and Mh+ = 182 GeV. The relevant decay branching ratios are 
B{A hh) = 0.90, B{H bb) = 0.90, B{H+ t+u) = 0.90, and 5(r+ ti+u) = 0.11. 
The relevant total decay widths are Ta = 5.6 GeV, Th = 5.5 GeV and Th± = 0.68 GeV, 
respectively. They are summarized in Table IH] for a quick reference. In this case, is 
about the same as Ma, and sin(/5 — a) ~ 1, hence, the production rates of AH~^ and HH^ 
are about the same. Since the decay branching ratios oi A ^ bb and H ^ bb are also similar, 
we shall include both the AH^ and HH^ event rates when we compute the signal rate in 
this case. Following the same analysis procedure as that for the small Ma case, we show the 
numbers of signal and background events in Table IVI After imposing all the kinematic cuts, 
the signal-to-background ratio increases by two orders of magnitude at the cost of half of 
signal events, and there are about 140 signal events and 50 background events. In Fig. lUJ we 
show the distributions of the invariant mass m^^ (after imposing the basic cuts, I^t and 
cuts) and transverse mass uit (after imposing all the kinematic cuts), which are obtained 
after smearing the observable parton momenta by 50%V^ to mimic the effect of the finite 
resolution of the detector. 

C. Ma = 250 GeV 

When the mass of the is larger than the sum of the masses of top and bottom quarks 
but not too large that the supersymmetric decay modes become significant, the dominant 
decay mode of is likely to be —>■ tb. Here, we study how to detect the signal event 
in this decay mode. 

To simplify the discussion, we shall concentrate on the semi-leptonic decay mode of top 
quark, and the bb decay mode of A or H. Hence, the signal events considered here are 
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TABLE V: Numbers of AH'^ and HH'^ signal and background events for Case B, with Ma = 
166 GeV and Mh = 163 GeV in the bbn'^^T channel, at the LHC with an integrated luminosity 
of 100 fb~^. The 6-tagging efficiency (50%, for tagging both b and b jets) is included, and the 
kinematic cuts listed in each column are applied sequentially. 





Basic Cuts 


> 50 GeV 


p'^ > 40 GeV 


155 < M^^ < 175 GeV 


With smearing 


AH+ 


126 


111 


85 


71 


65 


HH+ 


129 


114 


87 


72 


66 


Wbb 


11560 


3102 


840 


33 


34 


tb 


1221 


607 


164 


10 


10 


Wg 


783 


318 


11 


5 


6 


tt 


108 


79 


18 


1 


1 


Signal [S) 


255 


225 


171 


143 


131 


Bckg {B) 


13672 


4106 


1031 


48 


50 


S/B 


0.02 


0.05 


0.17 


3.0 


2.6 


s/Vb 


2.18 


3.51 


5.33 


20.8 


18.6 


+ B/S 


0.46 


0.29 


0.20 


0.09 


0.10 



produced via 

qq' A{^bb)H+{^bt{^be+ ly)), 

qq' bb)H+{^ bt{^ b £+ i/)) , (19) 

where = e+ or fi^. The signature of the signal event is four 6-jets plus one isolated lepton 
and missing transverse energy. To detect the signal event, we require the four 6-jets and 
the isolated charged lepton to be well separated in AR with large pr in the central rapidity 
region, in addition to a signature. In the following, we assume that the efficiency for 
tagging all four 6-jets is 25%, which is included in the calculation of the event rates to be 
discussed below. The dominant SM background for such kind of signature comes from the 
tibb events, produced via 

qq,gg -^bbt{^be+ iy)t{^bW-). (20) 
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FIG. 9: Distributions of invariant mass m^^^ (after imposing the basic cuts, and p!^ cuts) and 
transverse mass rriT (after imposing all the kinematic cuts) for Case B, with M4 = 166 GeV and 
Mh = 163 GeV, after smearing the observable parton momenta by 50% \/^. 



When the decay products of W~ escape detection, this background event will mimic the 
signal event signature. At the LHC, the tihb rate is dominated by the gluon fusion process 
due to its large parton luminosity. As shown in Table IVIl after imposing the following basic 
kinematic cuts:^ 

> 15GeV, |r/^^|<3.0, ^T>20GeV, 
p^T > 20 GeV, \ri^\ < 2.0, 
AR > 0.4, (21) 

the gg tihb rate is larger than the qq — > tihb rate by two orders of magnitude. Here, 
we have required the separation in Ai? between any two final state partons among the four 
6-jets and the charged lepton to be larger than 0.4. 

The signal rate depends on SUSY parameters. We take Ma = 250 GeV, tan (3 = 40, 
and all the other SUSY parameters (soft-breaking masses and fi parameter) to be 500 GeV. 
Consequently, Mh = 248 GeV and Mh+ =261 GeV. The relevant decay branching ratios 
are B{A bb) = 0.89, B{H bb) = 0.89, B{H+ tb) = 0.79, and B{t b e+u) = 0.11. 
The relevant total decay widths are Ta = 7.9 GeV, Th = 7.8 GeV and Th± = 4.2 GeV, as 



Here, we impose a set of stronger cuts than those given in Eq. H17|l in order to have a more reliable 
estimate of the ttbb background rate. 
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TABLE VI: Numbers of AH^ and HH^ signal and background events for Case C, with Ma_ = 
250 GeV and Mh = 248 GeV in the bbbbi'^^T channel, at the LHC with an integrated luminosity 
of 100 fb~^. The 6-tagging efficiency (25%, for tagging all the four b jets) is included, and the 
kinematic cuts listed in each column are applied sequentially. 





Basic Cuts 


Veto Cuts 


AH+ 


18 


18 


HH+ 


18 


18 


gg ttbb 


1632 


3 


qq -^■ ttbb 


105 





Signal (S) 


36 


36 


Bckg (B) 


1737 


3 



summarized in Table |n] For this set of parameters, the event rates of the AH~^ and HH^ 
modes are about the same due to the decouphng hmit, i.e., becomes large. Furthermore, 
the total signal rate (the sum of the AH^ and HH'^ rate) is about one hundredth of the 
tthh background rate, cf. Table IVH after the basic cuts, and the sum of the AH^ and HH^ 
signal rates yields about 36 events at the LHC. We have checked that our calculation of 
the tihb cross section agrees with that in Refs. 
evaluating parton distribution functions and strong coupling to be the invariant mass of the 
final state particles. (The ttbb cross section can vary as much as 50% with a different choice 
of the scale js^.) 

As noted above, the tihh background events contain an extra W~ boson in the final state, 
as compared to the signal event. To suppress this large background event, we need to veto 
the additional W~ boson which can decay into di-jets or a charged lepton plus neutrino. 
When W~ decays into the di-jet mode, there will be additional hadronic activities in the 
central rapidity region, hence, one needs to apply kinematic cuts to veto the additional 
high pt jets in the central rapidity region. When W~ decays into a charged lepton plus 
neutrino, an additional charged lepton is produced in the final state. Below, we first discuss 
the case that W~ decays into the leptonic mode with the charged lepton {^~) being e~ 




where we have chosen the scale for 
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or ij,", which accounts for two ninth of the top quark decay branching ratio. We veto the 
additional charged lepton with a transverse momentum exceeding 10 GeV and rapidity (in 
magnitude) less than 3.0. It turns out that the veto efficiency is so high (around 99%) 
that the background rate is down by three orders of magnitude and yields 2 events of 
t — >• bW~{^ z/£~) with fr = e~ or Consequently, the signal rate is about one order 
of magnitude larger than this background rate. Though the above analysis is adequate for 
£~ — e~ or it is less accurate for l~ — r" when t~ decays further to either leptons 
or hadrons. (The decay branching ratio of W~ vt" is about one ninth.) Hence, we 
improve the calculation of the t — > hW~{-^ Pt~) background rate by considering the leptonic 
decay mode t~ £~z/^z/ (with i" = e~ or ji^) and the hadronic decay modes of tau, 
such as — > Ti^Vr or p~Vt- For the £~z/^z/ decay mode, we veto events in which 

pfn > 10 GeV and \r]^ \ < 3.0. After taking into account the decay branching ratio of tau 
into the leptonic mode (about 35%), this amounts to 1 event of f ^ bW~{-^ Vt~{-^ £~Utu)) 
passing through the vetoing cut. It is straightforward to find out the background rate of 
t — > bW~{-^ Pt~(— > 7r~Ur)), after vetoing the additional central jet (i.e., 7r~ jet) that has 
transverse momentum larger than 10 GeV and rapidity (in magnitude) less than 3.5. We 
find that the vetoing efficiency for the t~ — » 7r~z/^ event, whose decay branching ratio is 
about 11%, is so high that its contribution to the tibb background rate is negligibly small. 
The similar conclusion also holds for the t~ — > UrP^i-^ tt^tt*^) mode, whose decay branching 
ratio is about 25%. For this decay mode, we first check whether the separation of tt" and 7r° 
in AR is within 0.4. If yes, we sum up their momenta to form one big jet and follow the same 
analysis as that for t~ TT'i/r- If not, we veto any additional central jet (i.e., vr" or 7r° jet) 
that has transverse momentum larger than 10 GeV and rapidity (in magnitude) within 3.5. 
In order for the t bW~{^ qq') background event to mimic the signal event, the decay 
products of W~ have to escape detection. When the separation of q and q' in Ai? is within 
0.4, we sum up their momenta to form one big jet and veto the event when the transverse 
momentum of the big jet is larger than 10 GeV and rapidity (in magnitude) smaller than 
3.5. Otherwise, we veto any additional central jet (i.e., q or q' jet) that has transverse 
momentum larger than 10 GeV and rapidity (in magnitude) within 3.5. This amounts to a 
negligibly small number of t — > bW~{-^ qq') events at the LHC, for its vetoing efficiency 
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is close to 100%. In Table IVIl we sum up all the decay modes of W~ discussed above 
to estimate the background rate ttbb bbbbi^I^T which is about one order of magnitude 
smaller than the signal rate. Furthermore, the distribution of the transverse mass of 
and ^T, which is the transverse mass of the boson in the signal event, can be used 
to further discriminate the signal from the background events. This is because the vetoing 
cuts generate additional missing energy in the background event, while in the signal event, 
the is entirely generated by the neutrino from boson decay. It is desirable to check 
the above estimate using a full event generator (including effects of parton showering and 
hadronization) combined with a realistic detector simulation, which is however beyond the 
scope of this paper. 

Next, we illustrate how to reconstruct the mass of A and in the signal event v4(— »• 
bb)H^{—^ bt{-^ b u)) . A similar method also applies to the HH^ signal event. Since both 

and v come from an on-shell boson decay, we can use the W-hoson mass constraint 
and f^T information to specify the longitudinal momentum {pz{v)) of the neutrino. The 
above procedure leads to two possible solutions of Pzi.^)- Sometimes, only one of the two 
solutions is kinematically allowed, but most of the time both of them are physical solutions 
for a signal event. Therefore, one has to fix a prescription to choose the one which 
will most likely give the correct distribution of the invariant mass of e"*", z/g and b. Since 
l^-boson, top quark and boson are all heavy particles, they and their decay products 
are likely to be produced in the central rapidity region. Hence, we choose the one with the 
smaller magnitude in Pz{^) to reconstruct the boson. In the signal event, there are 
four 6-jets. (We consider the case that a &-jet cannot be experimentally distinguished from 
a 6-jet.) We pair the reconstructed boson with any one of the four 6-jets to calculate 
the invariant mass of W^b and choose the 6-jet that yields a mass closest to the nominal 
top mass (174 GeV) to be the one produced from the decay of t. In Fig. [TUT a). we show 
the distribution of the reconstructed top quark mass which peaks around 174 GeV with a 
17 GeV half-of-maximum width. In order to pair the two correct 6-jets, among the remaining 
three 6-jets, to reconstruct the mass of the CP-odd Higgs boson A, we make use of the mass 
inequality implied by the MSSM mass relation, Eq. (0). Namely, the mass of A is smaller 
than the mass of . (As becomes large, the mass of H is about the same as the mass 
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FIG. 10: The invariant mass distributions of the 
boson , (c) CP-odd Higgs Boson A, and (d) 
Mh+ = 261 GeV, Ma = 250 GeV and Mr = 
discussed in the text. 
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of A.) We loop over the remaining three 6-jets (labelled as &i,2,3); and calculate the invariant 
mass of thi, denoted as mtb^, and the invariant mass of 62^3, denoted We choose 

the assignment of the three 6-jets so that rritb^ > rrib^b^ and the difference in m^^,^ and ^^263 
closest to M^. Following this procedure, we obtain the distributions of and Mh+ which 
are shown in Figs. ITUT b) and lTUT c). The distribution of the reconstructed Ma peaks around 
250 GeV with a 10 GeV, half-of-maximum width, and the reconstructed Mh+ peaks around 
261 GeV with a 20 GeV half-of-maximum width. The broad width in the reconstructed Mfj+ 
distribution is caused by the wrong assignment of the 6-jets, for the difference in Mfj+ and 
Ma (about 13 GeV) is at the same order as the decay widths of A and H^. Hence, even with 
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the crude analysis discussed above, it is possible to reconstruct the mass of A and in 
the AH~^ bbbt signal event. Similarly, Mh can also be reconstructed for the HH^ bbbt 
signal event, whose result is shown in Fig. [TUT d). A more complicated analysis using the 
likelihood method, such as the neural network analysis, would certainly improve the above 
conclusion. 

Finally, we comment on the effect due to the finite resolution of the detector. Again, we 
smear the four-momenta of the final state partons (except neutrino) by assuming a Gaussian 
error of 50%\/E, where E is the energy of the parton. Repeating the above analysis, we 
obtain the distributions of the reconstructed Mt, M4 (or Mh) and Mh+ as shown in Fig. ITTl 
Except the small difference in the height of the peak, they are similar to those obtained for 
a perfect detector, as shown in Fig. E3 

D. Discussions 

Before concluding this section, we would like to summarize the results of our Monte 
Carlo analysis and discuss the need for a realistic simulation of the detector and a full event 
generator analysis. In this work, we considered three different signal scenarios as given in 
Table ITTl In Case A and Case B, the dominant decay mode of is t^v^. For Case 

A, the SUSY parameters of the Higgs sector are in the mixing regime [sin(/5 — a) ~ 0.3-0.8], 
so that Mh differs from Ma- For Case B, the SUSY parameters of the Higgs sector are in 
the decoupling regime [sin(/3 — a) ~ 1], so that Mh is almost the same as M^. In Case C, 
the tb mode opens and becomes the dominant decay channel of , so that it provides a 
different event signature. The results of our Monte Carlo analysis are separately summarized 
as follows. 

To detect the signal events for Case A and Case B, we consider the decay modes oi A ^bb 
(or H — s> bb) and — > t^Vj. with r+ — n^Vj. in the AH^ (or HH^) events. The relevant 
decay widths and branching ratios are summarized in Table E] for a quick reference. The 
dominant SM background events for these cases come from Wbb, tb, Wg, and ti events, as 
defined in Eq. (fTBj) . 

Firstly, we discuss Case A, in which the production rate of HH^ is smaller than AH^ 
by a factor of sin^(/5 — a) ~ 0.11. After imposing the basic cuts, specified in Eq. (fTTjl . 
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FIG. 11: The invariant mass distributions of the reconstructed (a) top quark t, (b) charged Higgs 
boson H^, (c) CP-odd Higgs Boson A, and (d) CP-even heavy Higgs boson H, for Case C, with 
Mh+ = 261 GeV, Ma = 250 GeV and Mh = 248 GeV, after smearing the observable parton 
momenta by bOYoVE and imposing all the kinematic cuts discussed in the text. 



and assuming a 50% of 6-tagging efficiency for tagging both 6-jets in the event, we obtain 
the signal and background event rates as shown in the second column of Table IIIII for the 
LHC with an integrated luminosity of 100 fb^^. By noting that the background events 
(particularly the Wbb events) tend to have smaller missing energy, we impose a cut, > 
50 GeV, to suppress various background rates. The results are shown in the third column of 
Table III. Furthermore, the polarization of r"*" in the signal events (left-handedly polarized) is 
different from that in the SM background events (right-handedly polarized). Consequently, 
the typical transverse momentum of vr^ from the decay of r"*" is higher in the signal events 
than in the SM background events, cf. Fig. HI Hence, we impose another cut, p'!^ > 40 GeV, 
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to further suppress the background rates and increase the signal-to-background ratio. This 
result is shown in the fourth column of Table IIIIl To detect the AH^ signal event, we 
select the invariant mass of the bb pair so that it is within 10 GeV around M^. The result 
is shown in the fifth column of Table lllll As shown, after these cuts, the significance of 
the signal event increases by a factor of 5 at the cost of about half of the number of signal 
events. A similar result for detecting the HH^ signal event is given in Table IIVI Up 
to now, the analysis was carried out at the parton level with a perfect measurement of the 
observable parton momenta. In reality, the detection efficiency of the detector cannot be one 
hundred percent. To estimate the effect due to the finite detection efficiency of the detector, 
we smear the four-momenta of the final state observable partons (excluding neutrino) by 
assuming a Gaussian error of 50%V^, where E is the energy of the parton. After repeating 
the above analysis with the basic cuts, the ^t, Vti rrifji cuts, we obtain the signal 
event rates as shown in the last column of Table IIIIl and Table IIVI for the AH^ and HH^ 
channels, respectively. It is clear that the AH^ signal is much easier than the HH^ signal 
to be identified for Case A. Furthermore, from the distribution of the hb invariant mass 
771^5, one can extract out the A boson mass Ma from its peak position, cf. Fig. From 
the distribution of the transverse mass of defined as Eq. (jl8|) . one can extract out the 
H'^ boson mass Mh+ from its jacobian peak position, cf. Fig. IHl Given the measured Ma 
and Mii+ values, one can then test the mass relation, Eq. (^, to confirm the MSSM Higgs 
sector. This is the conclusion for the case that a AH^ signal, as predicted by Case A, is 
found. If it is not found, then one can readily constrain the MSSM. For example, according 
to the last column in Table IIIIl if a AH^ signal is not found, then at the two standard 
deviation level, the product of B{A bb) and B{H^ J^tT^) in the MSSM has to be less 
than X (0.91)(0.98) = 0.082 for Ma = 101 GeV.^ We note that, as discussed in the 

previous sections, this upper bound on B{A —>■ bb) x B{H^ ^t't^) as a function of Ma 
(i.e., rriiji) is not sensitive to the other SUSY parameters in the MSSM. 

Surely, the above bound can be improved by including the negatively charged mode 
AH~ of the signal event (whose rate is about half of the AH^ rate, cf. Fig. and the other 

^ Here, we require the number of signal events to be less than two times the statistical fluctuation {^/~B) of 
the SM background events. 
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decay mode of t"^ (such as r"*" p^i^r and r~ —>■ p~Vt)- Since the decay branching ratio 

B(r+ — s> P^Vt) is a factor of 2 larger than B(r+ — vr+z/^), and the dominant decay mode 

of p+ ^ tt+tt'^ produces a tt^ whose typical momentum is smaller than that produced in 

r+ n^Ur, we expect that including the r+ p^v^ channel would increase the event rate 

by about a factor of 2. Therefore, after including both the AH^ and AH~ channels with r 

decays into the vrz/ or pv modes, the signal rate will increase by about a factor of 3. Assuming 

the background rates increase by the same factor, the resolution power will increase by a 

factor about v^- For example, the above constraint will read as B{A bb) x B{H^ — »• 

z/^r±) < 0.08/^3 = 0.046 for Ma = 101 GeV in the MSSM. However, a reliable conclusion 

can only be drawn from the study using a full event generator (that predicts the distributions 

of final state hadrons) folded with realistic detector simulation. 

Secondly, we consider Case B, in which the production rates of HH^ and AH^ are about 

the same. Following the same analysis procedures as in Case A, we obtain the results shown 

in Table |3 and Fig. IHl According to the last column in Table if a AH^ signal is not 

found, then at the two standard deviation level, B{A —>■ bb) x B[H^ i^tT^) in the MSSM 
2^/50 

has to be less than — — x (0.90) (0.90) = 0.18 for Ma = 166 GeV. After including the AH' 
65 

production channel and t ^ pv decay modes, we expect the above bound to improve by 
about a factor of ^3, i.e., B{A bb) x B{H^ u^r^) < 0.18/^3 = 0.1 for Ma = 166 GeV 
in the MSSM. 

Thirdly, we consider Case C, in which predominantly decays into the tb mode. For 
the CP-odd Higgs boson decays into a bb pair, the dominant SM background comes from the 
bbti production when the decay products of , produced from i decay, escape detection. 
The sensitivity of the LHC to this case is presented in Table IVll and Figs. ITUlandlTTl Due 
to the complexity of the background simulation, the result of our parton level Monte Carlo 
analysis for the background event can only be viewed as an estimate. A more reliable 
calculation using the full event generator with detector simulation is needed to draw a 
definite conclusion. Nevertheless, we illustrated that the tb decay mode of the signal event 
can be detected at the LHC with a reasonable resolution on the determination of Ma and 
Mfj+ , cf . Figs. ^1 and ^2 A crude estimate from Table IVII reveals that if a AH^ signal is 
not found, then at the two standard deviation level, B{A ^ bb) x B{H^ —>■ tb) has to be 
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TABLE VII: Numbers of AH^ signal and background events for Case A, with M4 = 101 GeV 
in the bbir^l^T channel, at the Tevatron with an integrated luminosity of 30fb~^. The 5-tagging 
efficiency (50%, for tagging both b and b jets) is included, and the kinematic cuts listed in each 
column are applied sequentially. 



AH+ 
HH+ 
Wbb 

tb 

ti 

Wg 



Basic Cuts 



11 
1 

383 
6 





> 50 GeV 



27 
3 





p J > 40 GeV 



90 < Ma < 110 GeV 
or 

105 <Mh < 125 GeV 
3 








less than ^— x (0.89) (0.79) = 0.135 for Ma = 250 GeV in the MSSM. 
18 

Finally, we comment on the potential of the Fermilab Tevatron. As shown in Fig. ^ the 
production rate of the AH~^ event at the Tevatron is only sizable for small M^. Assuming 
an integrated luminosity of 30 fb~^ at the Tevatron, the event yield for Case A is shown in 
Table IVllI It is evident that it is challenging to detect such a signal at the Tevatron. 



IV. CONCLUSIONS 



Most production processes studied in the literature for testing the MSSM contain at 
least two SUSY parameters (such as tan j3 and Ma) in the search for supersymmetric Higgs 
bosons. Furthermore, the detection efficiency of the signal event depends on the assumed 
decay channels of the SUSY particles, hence, on the detailed choice of SUSY parameters. 
If the signal is not found after comparing experimental data with theory prediction, it is 
a common practice to constrain the product of the production cross section and the decay 
branching ratios of final state SUSY particles as a function of the multiple-dimension SUSY 
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O Case A (X = bb, Y = T*v) 
□ Case B (X = bb, Y = l\) 
V Case C (X = bb, Y = tb ) 
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FIG. 12: Constraints on the product of branching ratios B{A bb) x B{H^ t^i^t) as a function 
of Ma for Case A and Case B, and B{A bb) x B{H^ tb) for Case C, at the LHC, where t+ 
decays into tt^i't- channel. 

parameter space of the MSSM. 

In Ref. j2^, a novel proposal was made to study the AH"^ production process at hadron 
colliders. It was pointed out that this process possesses the following interesting properties: 

• Its Born level production rate depends only on one SUSY parameter that can be 
determined by kinematic variables (e.g., the invariant mass of the bb pair). 

• Its higher order production rate is not sensitive to detailed SUSY parameters through 
electroweak radiative corrections. 

• Its final state particle kinematics can be properly modelled without specifying any 
SUSY parameters. Hence, the detection efficiency of the signal event can be accurately 
determined. 

• If the signal is found, it can be used to distinguish the MSSM Higgs sector from its 
alternatives, e.g., THDM. 

• If the signal is not found, one can constrain the MSSM by limiting the product of 
decay branching ratios alone, without convoluting with the production cross section. 



1^ 
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At the LHC, the AH~^ signal event can be produced from the tree level process qq' 
W^* AH^, whose rate depends only on Ma and Mh+, for the coupling of W-A-H^ 
is fixed by the weak gauge coupling g, required by SU{2) gauge invariance. Furthermore, 
in the MSSM, Ma and Mh+ are related as Eq. therefore the AH^ production rate 
depends only on one SUSY parameter, i.e., M4, which can be kinematically determined 
from experimental data by reconstructing the invariant mass of the decay particles of A. 
For example, when A decays into the bb channel, the invariant mass of bb pair reveals the 
mass of A. Since both A and are scalar particles, the angular distributions of their 
decay particles can be accurately modelled, i.e., an isotropic distribution in its rest frame. If 
the signal is found, the mass relation, Eq. (^, can be tested to distinguish the MSSM from 
the THDM in which this mass relation does not generally hold. If the signal event is not 
found, then one can constrain the MSSM parameters by limiting the product of the decay 
branching ratios B{A ^ X) x B{H^ ^ F) as a function of Ma, where X and Y are the 
decay channels of A and H^, respectively, studied by experimentalists. 

In Sec. II, we show that the electroweak radiative corrections that depend on detailed 
SUSY parameters are usually small in comparison to the uncertainty in higher order (beyond 
the NLO) QCD corrections, parton distribution functions, or the accuracy of the experimen- 
tal measurement. 

To test whether such a signal can be detected at the LHC, we performed a Monte Carlo 
study at the parton level in Sec. III. We consider three cases of MSSM to cover the decay 
modes of A{—>- bb)H^{-^ t^i^t) and A{—>- bb)H^{-^ tb). We concluded in the last part of 
Sec. Ill that at the LHC this signal event can indeed provide useful information about the 
MSSM Higgs sector. For example, if the AH^ signal is not found in the decay mode of r+ — >■ 
vr+z/^, then we can constrain the product of branching ratios B{A bb) x B{H^ t^i't) 
as a function of Ma-, as shown in Fig. ^| This corresponds to Case A or Case B defined 
in Table |H] In case C, for Mh+ > rrit + nif,, not finding the signal event implies an upper 
bound on B{A bb) x B{H^ tb) for a given Ma-, and therefore constrains the MSSM as 
a function of one SUSY parameter. Including the negatively charged channel AH~ and the 
pv decay mode of r can tighten the above bounds roughly by a factor of -\/3- However, to 
have a more accurate conclusion, a full event generator with detector simulation should be 
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FIG. 13: The LO (dotted lines) and NLO QCD (solid lines) cross sections of the AH+ and AH' 
pairs as a function of at the VLHC (a 200 TeV pp cohider) . 

used to repeat the analysis outlined in this paper. 

From Fig. ^ we see that the AH~^ rate becomes very small (less than about 0.1 fb) at 
the LHC once Ma_ is larger than 400 GeV. Hence, to cover the whole mass spectrum of the 
TeV scale MSSM, we should find a high energy collider that can be sensitive to this process 
for Ma approaching the TeV region. This could be one of the motivations for proposing a 
future Very Large Hadron Collider (VLHC), a 200 TeV proton-proton collider. The signal 
rates at the VLHC are shown in Fig. [T^ 
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V. APPENDICES 

A. Renormalization 

In this paper, we adopt the on-shell renormahzation scheme developed by Dabelsteinj2^ 
to calculate the one-loop electroweak corrections. The standard model parameters are fixed 
by defining aem, i^w and m^, and the additional SUSY parameters in the Higgs sector^ are 
fixed by the following renormalization conditions: (1) the tadpole contributions {Th^ = 0, 
= 0)) (2) the on-shell condition for the mass of A, (3) the on-shell condition for the 
wavefunction of A, (4) a renormalization condition on tan (3 (which requires 5vi/vi = SV2/V2), 
and (5) a vanishing A — Z mixing for an on-shell A. 

B. Calculation of F(i)(g2) 

The one-loop correction to the renormalized form factor of the W^H'^A vertex, apart 
from the effective weak gauge coupling g, can be written as 



= Z]/^Zljl^_{l + 6F^HA 

+F^lAMlMl,q^)]-l, (22) 

where Zaa and Zh+h~ are the finite wavefunction factors for the renormalization of the 
external Higgs bosons A and H^. In our scheme, 

Zaa = 1, (23) 
Zh^h- = 1 - n'^+^_ (Ml + m^) + U'^AM'a) , (24) 

where n^^(MD denotes taking the derivative of the two point function n^^(A;^) of the 
CP-odd scalar A with respect to k"^ at k'^ = M\, etc. The terms inside the curly bracket of 
Eq. ()22|1 arise from the renormalized vertex function of WHA. Fy^jj j^{p\, p]j ^ q^) represents 

^ There are 7 parameters in the Higgs sector of the MSSM. They are 17', vi, V2, mi, 7712, and 7713. Beyond 
the Born level, the wavefunction renormalization factors Z^^ and Z^^ also need to be introduced to 
renormalize the theory, where Hi and H2 denote the two Higgs doublets in the model. 
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the one-loop contribution of the one-particle-irreducible (IPI) diagrams with p^pf^jq"^ as 
the four-momentum square of the incoming A, and particles, respectively. 5Fwha 
is the counterterm contribution resulting from the field renormalization of and A: 

H^A ^ H^A (^1 + ^5Zh+ + ^5Za^ ■ (25) 

In terms of the independent counterterms fixed by the renormalization scheme, the wave- 
function counterterms SZh+ and 6Za can be written as (sin^ + {cos^ P)6Zh2 which 
is found to be equal to — |n^y^(M^). We note that in 6Fwha the contributions from the 
counterterms of the weak gauge coupling and the wavefunction renormalization of the W- 
boson are not included, because they should be combined with the VT-boson self energy 
contribution to derive the running weak gauge coupling g{q^). In our numerical calculation, 
we use 

f = ^^ml^Gp . 

In summary, the one-loop electroweak correction to F^^\q'^) is found to be 

-^n'^+^_(Ml + ml,) - \ll'^AMl). (26) 
In the above equation, the top- and bottom-loop contribution to F^jj^ is given by 

F'^ir'-'iq^pWn) 

= E ^^'Li<i'yAyH), (27) 
fff'=ttb,bbt 

with 

^whakPa^PhiQ. ) — ^iQ^2 yf yPA^^zi Ph'^32 

+ {2pA ■ PH - pI)C!J^' - {2PA ■ PH - pI)CU^' 

+{D + 2){Cii^' - CiP') + p\Ci(^' - {2pA ■ PH 

+pl)Ci(^' - 2p\Cii^' -{D- 2)ClF - ^)CW' 

- {q' + m)) C//^'} - cs-^^VsVrmfms.Cl^r (28) 

where c/ = +1 and —1 for ///' = tth and respectively, and are defined in terms 

of the Passarino-Veltman functions [s^ with 

CiP' = Cij {p%Ph, (pa + Puf] ruf, mf,mf) . (29) 
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The stop- and sbottom-loop contribution is given by 

N 2 

plPI(squark) / 2 2 „2x _ /o 
^WHA KPaiPhi^ ) — lQj^2^ / > 

j,j,A;=l 

-VIj^Du bj, A]) X[b*, ik, H-]C'^'^'^ } , (30) 

where Un, Dn arc the rotation matrices for stops and sbottoms between the weak eigenstate 
basis (/ — L, R) and the mass eigenstate basis (i = 1, 2), respectively. A[/*, /j, 0^] represents 
the coefficient of the /f/j^fe interaction in the MSSM Lagrangian, as hsted in Appendix C, 
and 

cM^f^ = (Cn-Cu) {pl,pl,q';mf^,mf,,mj.) . (31) 

The quark (top and bottom) and squark (stop and sbottom) loop contributions to the self- 
energies IIaa{(i'^) and IIh+h-{q^) can be found in Appendix E. 



C. Squark couplings with and A 

The mass eigenstates of the squarks are related to their weak eigenstates by the rotation 
matrix O-j with /j = Yli Opifi^ where i — 1,2 and I — L,R; Oj- — Un and Dn for f — t 
and b, respectively. In terms of the mixing angles 9f, we have 

f fA J cose, .n0,\ f h\ p,, 
\h J \-sme, cosO/ J \ fK J 

Here, we define the mixing angle ^/ so that /i is hghter than /2. 

The coupling constants among the weak-eigenstate squarks and the Higgs bosons are 
defined through the Lagrangian 

>c = • • • + A[/;, fj, 0, ...]/;/>, ... + •• •. (33) 

Hence, the coupling constants for the mass-eigenstate squarks are a linear combination of 
the couplings for the weak-eigenstate squarks, and 

mjl<p,-] = mo{^,o{;f'j,4>,...] 

= (34) 
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The relevant couplings A[//, f'j, 0, ...], denoted as X^*^,^^ ., are hsted below. 

-\/2 

HitLH- ~ ("^VK 2/3 — ml tan /3 — m^ cot /3), 

. \/2mt , . „ . 

hliaH- = (A cot /3 + //) , 



Am 



V 

2\/2mtmb 



V sin 2/3 ' 
^flhA = ^rjRA = 0, (/ = 6), 

\bRA = (^6tan/3 + /i), 

. imt , . „ . 
At* IrA = [A cot /3 + /x) , 

2 2 

^bibi^AA = ^ tan^ /3 + ^(T, - Q^s^) cos 2/3, 

777^ (7^ 

hll^AA = cot^ /3 + ^(Tt - Qtsl,) cos 2/3, 

777^ 

\bRAA = —jr cot^ ^ + x^*^^ 

hlbLH^H- = tan^ /3 + ^(T, - g,s^) cos 2/3, 

2 



h*~bRH+H- = — tan /3 + — QbSvy cos 2/3, 

^ V Zi 



_2777^ (7^ 

At*f^/^+,/- = ^^tan2/3 + ^(T,-Qts^)cos2/3, 
MiIrH+h- = cot^ /3 + ^Qts^ cos 2/3, 

where T^, T^, and are |, | and respectively, and 

\*bjH+ = \*tjH-^ 

for I , J — L, R and f = i,b. 
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D. Squark contributions to the p parameter 



The squark one-loop contribution to the p parameter is given by 

Ap = p - 1 = -4v^Gi.Re[An^i(0) - Anp(O)], (52) 



with 



f=t,bi,j=l 



xB{q'^;m^j^,my), (53) 

N ^ 

"^"^'^^'^ = 3^ 5: \Uu\'\Di/Biq^;ml,ml), (54) 

where O^^ = Uu and 0'}^ = Dn; B{q^;ml,ml) = A^mj) + A{ml) - 4522(g^ , m|). By- 
using the expression 

B{0;ml,ml) = --{mj + ml) + 7^"^% In ^ . (55) 

Eq. (0) is deduced under the assumption that = Mq ~ ^ M|, ^ and the 
stop mixing is large (m(|Xt| ~ and m(,|Xfe| ~ 0), so that m^^ ~ 0{mz), which yields 
rrif ~ ^/2M, and rrir ~ mr ~ M. 



E. Self energies 

The (top and bottom) quark- loop contributions to the self-energies 11^^(0'^) and 
IIh+h-{(1^) in -D-dimensions are expressed in terms of the Passarino-Veltman functions jSTj 
by 

f=t,b 

+52l(g^ m/,m/)) + DB22{q'^ ,mf,mf) 

+mJ5o(g^m/,m/)} , (56) 



+B2i{q'^,mb,mt)) + DB22{q'^,mb,mt)} 



■Y^4:ybytmbmtBo{q'^,mb,mt). (57) 
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The stop- and sbottom-loop contributions are given by 

f=i,b «.i=i 

x\[f:J„A]X[f;J,,A]Bo{q',mj^,mj^) 
N, 



2E (58) 



167r2 

f=i,b *=i 



2 



«J=1 



16,,2E (59) 

f=i;b i=i 

The self-energy n^j;(/(g^) of the W boson was aheady presented in the hterature. For exam- 
ple, the quark-loop contribution can be found in Ref. |38|, and the squark-loop contribution 
in Ref. m. 
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